Rough Surfaces: is the dark stuff just shadow? 


“Who knows what evil lurks in the hearts of men? The Shadow knows!” ! 


Jeffrey N. Cuzzi', Lindsey B. Chambers’, Amanda R. Hendrix? 
‘Ames Research Center, NASA; ?University of California, Santa Cruz; ?PSI, inc. 


Revised, November 7, 2016 


Abstract 


Remote observations of the surfaces of airless planetary objects are fundamental to inferring 
the physical structure and compositional makeup of the surface material. A number of forward 
models have been developed to reproduce the photometric behavior of these surfaces, based on 
specific, assumed structural properties such as macroscopic roughness and associated shadowing. 
Most work of this type is applied to geometric albedos, which are affected by complicated ef- 
fects near zero phase angle that represent only a tiny fraction of the net energy reflected by the 
object. Other applications include parameter fits to resolved portions of some planetary surface 
as viewed over a range of geometries. The spherical albedo of the entire object (when it can 
be determined) captures the net energy balance of the particle more robustly than the geomet- 
ric albedo. In most treatments involving spherical albedos, spherical albedos and particle phase 
functions are often treated as if they are independent, neglecting the effects of roughness. In this 
paper we take a different approach. We note that whatever function captures the phase angle 
dependence of the brightness of a realistic rough, shadowed, flat surface element relative to that 
of a smooth granular surface of the same material, it is manifested directly in both the integral 
phase function and the spherical albedo of the object. We suggest that, where broad phase angle 
coverage is possible, spherical albedos may be easily corrected for the effects of shadowing using 
observed (or assumed) phase functions, and then modeled more robustly using smooth-surface 
regolith radiative transfer models without further imposed (forward-modeled) shadowing correc- 
tions. Our approach attributes observed “powerlaw” phase functions of various slope (and “linear” 
ranges of magnitude-vs.-phase angle) to shadowing, as have others, and goes in to suggest that 
regolith-model-based inferences of composition based on shadow-uncorrected spherical albedos 
overestimate the amount of absorbing material contained in the regolith. 


'From the 1930s radio serial “The Shadow”, originated by Street and Smith Publications, inc. (Stedman 1977) 


1 Introduction 


Inference of surface composition from remote sensing is a core tool in planetary science. In studies 
of planetary rings, there is added motivation to obtain quantitative results for compositional fractions 
because the relative abundances of icy and nonicy material can be diagnostic of meteoroid bombard- 
ment and ring age. Classical treatments (Harris 1961) show how the phase function of an object plays 
a key role in relating its observed brightness at some arbitrary phase angle (generally called I/F’, 
as discussed below) back to either its geometric or “spherical” albedo, either of which can constrain 
composition (see next section for definitions). If the phase function is not measured (generally the 
case for remote objects not visited by spacecraft), then representative phase functions (for, e.g.., inner 
and outer solar system, dark and bright, objects) can serve as substitutes, with some uncertainty. Most 
compositional modeling of this kind has been applied to low phase angle observations where surface 
shadowing plays a small role (see section 1.5) and low reflectivity is primarily due to absorbing mate- 
rials mixed into the near-surface regolith (eg. Emery and Brown 2004, Cruikshank er al. 2005). These 
observations are transformed into geometric albedos, the reflectivity at zero phase angle, and model 
inferences are often not even corrected for shadowing. 

The need to incorporate roughness in compositional modeling is certainly not unknown to afi- 
cionados of the technique. Detailed algorithms for calculating these effects have been known for 
decades (Lumme and Bowell 1981a,b, Hapke 1993, his chapter 12). The former calculate the shad- 
owing for single scattering separately in porous surfaces and macroscopic “holes”, adding in a model 
for multiple scattering in both structures. The latter assume a faceted surface structure and calculate 
a shadowing correction which is applied to an independently modeled smooth-surface albedo, and 
neglect multiple scattering. Both are forward models, based on assumed parameters. For a broad dis- 
cussion and background of planetary surface photometry and the effects of roughness and shadowing, 
the reader is also referred to Shkuratov et al. (2005, 2012), Hapke (2013) and Shkuratov et al. (2013), 
as well as Mushkin and Gillespie (2006), Jehl et al. (2008), Muinonen et al. (2011), and Goguen 
et al. (2010) which cover a representative selection of approaches and targets. In a few recent studies, 
extensive grid-based, Bayesian, or genetic algorithms have been used to separate roughness and com- 
positional effects (Cord et al. 2003, 2005, Souchon et al. 2011, Pilorget et al. 2016, Domingue et al. 
2016). 

In the planetary rings literature, attention to the effects of shadowing on particle albedo has 
been mixed. Shadowing was untreated in spherical albedo modeling by Doyle et al. (1989), Cuzzi 
and Estrada (1998), Bradley et al. (2010, 2013), Elliott and Esposito (2011), and Poulet and Cuzzi 
(2002). Poulet et al. (2003) worked with the near-zero phase angle reflectivity of a ring particle, after 
removing optical depth and geometry effects, and also neglected roughness treatments; they needed 
to add some 5-10% carbon to the ring material to match the overall albedo. Poulet et al. (2002), also 
working over a small range of low phase angles, explicitly treated ring particle roughness using both 
the Hapke (1993) and Shkuratov et al. (1999) models, finding the A and B ring particles to be rougher, 
but did not go so far as to derive particle compositions. 

Ciarniello et al. (2011) and Filacchione et al. (2013) normalize spectra of icy moons in an at- 
tempt to remove the effect of roughness, and make use of the secondary effect of internal absorbing 
material on depth of ice absorption bands (ring spectra were not quantitatively analyzed for composi- 
tion), and Hedman et al. (2013) also work with moderate phase angle reflectivities but do not account 
for surface roughness. It is possible, but has not been demonstrated, that modeling spectral band 
shape and contrast alone, as in Ciarniello et al. (2011), Clark et al. (2012), Filacchione et al. (2013), 
and Hedman et al. (2013) can eliminate the effects of surface roughness, but given the number of pa- 
rameters available to absorb the effects, including continuum absorber and regolith grain properties, 


one remains concerned (see end of section 2). For example, Hedman et al. (2013) analyzed moderate 
phase angle brightness values using a Shkuratov model for a smooth, flat surface, inferring significant 
radial variations in “effective” grain size that are strongly correlated with local optical depth varia- 
tions across the rings. However, both on-surface shadowing and interparticle shadowing can affect 
reflectance at moderate phase angles (Déau 2015), and either or both could vary with local optical 
depth in a ring because of local collisional dynamics. As discussed in section 2, further modeling is 
needed to separate radial (optical-depth-related) variations in the “effective” regolith grain size (Hed- 
man ef al. 2013) from variations in the macroscopically rough nature of the particles themselves, 
and the ensuing shadowing effects, or even optical-depth-dependent divergence in layer scattering 
properties from classical treatments. 

The present paper shows, in a simple way, how shadowing effects on spherical albedo can be 
substantial, but can also be constrained by observational determination of the integrated particle phase 
function. Such a treatment could remove one level of ambiguity from subsequent analyses, without 
relying on uncertain forward models, using only knowledge of the particle phase function that comes 
out of the same analysis that leads to the spherical albedo (see section 1.4 for more discussion). The 
approach may also be useful in modeling composition of NEOs, where phase functions are well- 
determined, or even main belt asteroids in many cases (Shevchenko et al. 2016, Penttila et al. 2016). 

Overview of the paper: After some basic definitions we discuss the reflectance of locally flat, 
smooth surface elements, and define key parameters and functions (sections 1.1 and 1.2). In section 
1.3 we introduce the role of surface roughness and the associated shadowing, and define a key param- 
eter to describe it. In section 1.4 we discuss phase functions for three dimensional particles covered 
with locally “flat” but “rough” surfaces as discussed above, and show that commonly observed “‘pow- 
erlaw” phase functions are a natural outcome of our simple assumptions.” In section 1.5 we discuss 
the planar albedos of rough flat surfaces and the spherical and geometrical albedos of particles cov- 
ered by them, and show how to convert an observed spherical albedo, for a rough and shadowed 
object, to the spherical albedo of an object covered by a smooth, locally flat surface made of the same 
regolith material (which traditional regolith models can address with less ambiguity). In section 2 we 
compare our simple model to the full roughness-shadowing forward model of Hapke (1993), finding 
decent agreement with less dependence on initial assumptions, and mention possible implications for 
spectral modeling. 


1.1. Basic definitions 


The reflectance of a resolved surface is typically expressed as I/F’, which is the ratio of the observed 
intensity J (ergs cm~? str~'sec~') to that of a perfect Lambert disk illuminated at normal incidence 
by a flux density 7F (ergs cm~? sec~'). Cuzzi (1985) presented a discussion and table relating the 
various photometric properties of flat Lambert and Lommel-Seeliger (LS) surfaces, and the albedo 
and phase function of a sphere covered with such surfaces. We first define the planar flux albedo 
Apr, or ratio of the total flux reflected by a flat surface (integrated over emergence angles e) to 
incident flux at incidence angle i, as: 
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*In this paper we discuss macroscopic objects covered by regolith. For clarity we will refer to the macroscopic objects 
as particles and the components of their regolith as grains, and use different symbols for their respective albedos. 

3Hapke (1993) refers to App as Ap (hemispherical albedo), or plane albedo, or directional-hemispherical reflectance 
(his sections 10.E.4 and 10.D.2). 


Figure 1: Typical illumination and viewing geometry of a locally flat surface element with surface 
normal 7. Solar illumination is at incidence angle 7, and the viewer is at emission angle e, where 
[to = cost and ys = cose. The phase angle a is the angle between the incident and scattered ray, and 


is given by cosa = [lo + cosby/1 — p?s/1 = p?. 


where ju = cos(e), Wo = cos(z), and ¢ is an azimuth angle (figure 1). Note that App is not the same 
as geometric albedo p (Harris 1961), called A, by Hapke (1993), which is simply the //F at zero 
phase angle (Cuzzi 1985). We define the essential common underlying parameter of our reflectance 
functions as R, the planar flux albedo of a macroscopically smooth, flat, surface element illuminated 
at normal incidence. For instance, a Lambert surface with Apr = FR has an angle-independent 
reflected intensity [/F' = Ry. We chose this definition of R because it is independent of surface 
physical structure and can be measured in the laboratory. 

A second important albedo we will use extensively is the spherical albedo A, of an entire 
particle (also known as the single-scattering albedo): 
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For smooth surfaces (table 1), there is a simple relationship between 7 (the parameter most cleanly 
constraining composition) and A, (the parameter most cleanly characterizing particle brightness over- 
all). In general however, it is necessary to retain them both. 


1.2 Smooth flat surfaces 


Popular smooth, flat surface scattering laws include the Lommel-Seeliger (LS) law (Table 1) and the 
Lambert law, which is a subset of the Minnaert law: 
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where &; is the Minnaert parameter and a Lambert surface has & = 1. A Lambert law can be theoreti- 
cally demonstrated to be the high-reflectivity limit of rough particles; see van de Hulst (1957, section 
8.42), Lumme and Bowell (1981a), Shkuratov et al. (2005), and Muinonen et al. (2011). ALS surface 
has close connections to the single-scattering solution for a traditional diffuse or “classical” many- 
particle-thick layer, and it has been used often over the years for modeling the moon and various 
other dark objects (Hapke 1993, Veverka et al. 1989, Squyres and Veverka 1981). Counterintuitively, 
the planar flux albedo App(j1.) of a macroscopically smooth LS surface with high assumed normal 
reflectivity Apy(1) = R can exceed unity for small 1, (see table 1), and the single scattering albedo 
of a particle covered with such a surface is A, = 42/3, which can also exceed unity for R > 3/4 - 
both nonphysical implications. The fact that A, > R for a LS surface becomes more intuitive when 
one recalls that such a surface does not display its maximum reflectivity at normal viewing. More- 
over, since the LS law is the low-reflectivity, single-scattering limit for a layer of isotropic scatterers, 
clearly it would be inconsistent to assign a high reflectance R > 3/4 to a LS surface; thus nonphysical 
magnitudes of App and A, should not be encountered. While a highly reflective surface will simply 
not have the directional properties given by the LS expression, and is best modeled using a Lambert 
or Minnaert function, the smooth-surface behavior of dark objects should probably be based on the 
LS law instead of Minnaert surfaces with k < 1 for this reason (see section 1.5). 

In figure 2, we compare the Lommel-Seeliger law with two different scattering laws that are 
more compatible with bright surfaces. We plot J/F as a function of yu (the cosine of the emission 
angle) for three different values of 41. Dotted lines represent the Minnaert functions. Solid lines 
represent the classical multiple-scattering result for a semi-infinite layer composed of isotropically 
scattering monomers: I/F’ = S({1o, 1, @) /4p4 = R( po, UW, P)Mo/4(te + fo) (Chandrasekhar 1960; see 
figure 1). Using the notation of Hapke (1993), where w, is regolith grain albedo: 
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where H (Wee) is Chandrasekhar’s H-function, and for numerical calculations we will adopt the 
closed-form expression from Hapke (1993, equation 8.57) with x = yu or Uo: 


H(x) = ft - 20-9) r+ (1 are) in(442)]| (5) 
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where y = ,/1—w, andr, = (1 — y)/(1+7) (see section 3 for references to these parameters). 
Different solid curves in figure 2 are for regolith grain albedo w, = 0.3, 0.6, 0.9, 0.99, and 1.0. It 
can be seen that, for very low values of wz, the classical solution approaches the LS law (dashed line 
assuming R = 0.04). It can also be seen (especially for lower ju.) that the classical solutions are 
generally limb bright (brighter as 4 — 0), while the Minnaert solutions (dotted curves) are always 
limb dark; for the Lambert case (k = 1) the reflected brightness is independent of y. Minnaert 
solutions with & ~ 1.1 are in reasonable agreement with classical solutions for high grain albedos (left 
and middle panels); however, classical solutions for low incidence angles remain more limb-bright 
(right panel). In no case would a value of k > 1.15 or so be a reasonable match, overall. Minnaert 
surfaces for k < 1 share the properties of the Lommel-Seeliger surface, and are also nonphysical for 
highly reflective surfaces. 


1.3 Rough Flat surfaces: 


The above discussion applies to macroscopically smooth, unshadowed surfaces. Unfortunately, even 
for k > 1, Minnaert (or classical) layer reflectivities by themselves provide a poor match to the 
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Figure 2: Solid lines: Classical H-function reflectivities [/F for five values of regolith grain albedo 
Wy = 0.3, 0.6, 0.9, 0.99, and 1.0 (from bottom to top). Dotted lines: Smooth Minnaert surface re- 
flectivities for R = 1 and k = 1.0, 1.05, 1.15 (k = 1, the Lambert case, is constant). Dashed line: 
Lommel-Seeliger reflectivity for R = 0.04. 


variation of brightness with phase angle a (see figure 1) which have been observed for realistic solid 
objects (eg. Dones et al. 1993, Pitman et al. 2010) or laboratory surfaces (Helfenstein and Shepard 
2011, Shepard and Helfenstein 2011). 

The most obvious physical attribute lacking from Minnaert and classical layer models is signif- 
icant macroscopic relief on the target’s surface - facets, lumps and holes - and the shadowing thereby 
created. Relative to a smooth surface, an observer at higher phase angles sees more shadows, and 
thus a more sharply backward peaked phase function (section 1.4). Harris (1961) presented a version 
of the Minnaert law with a-dependent 7 and & in order to match lunar photometry quantitatively 
(see also Minnaert 1961). Shkuratov er al. (1999) and Kreslavsky et al. (2000) adopted a shadowing 
function of the form f(a) = exp(—Sa), originally suggested by Akimov (1980), as a prefactor to 
the smooth-surface reflectance laws, where S is a roughness parameter. Squyres and Veverka (1981) 
employed an empirical multicomponent f(a) factor, of which the dominant term is also of the form 
exp(—Sa). Hapke (1993) (chapter 12) models this effect from basic geometrical shadowing princi- 
ples (under various simplifying assumptions), obtaining functions of ju, and js multiplied by terms 
with more complicated functions of a. Goguen ef al. (2010) related S (in our notation) to a different 
physical model of slopes, and Shkuratov et al. (2005) models blocky lumps (for more discussion see 
these articles and references therein). None of the roughness/shadowing models can really distinguish 
the lengthscale over which roughness and shadowing acts, but comparisons of theoretical and exper- 
imental work (Shepard and Helfenstein 2011, Helfenstein and Shepard 2011) suggest the relevant 
scale is small - mm to cm. 

We find the heritage and simplicity of this form of prefactor shadowing correction appealing, 
and in fact below we show that, in combination with a flat facet reflectance law as simple as Lamber- 
tian, it captures much of the observed behavior of ring particle (and icy moon) phase functions. Here, 
we will explore two forms: 


f(a) =exp(—Sa) and f(a) = exp(—Sy/tan(a/2)). (6) 


We will refer to the first expression in equation 6 as the “simple exponential” form. The second 
expression is functionally similar to, but much simpler than, that derived by Hapke (1993), which 
matches the best-observed phase functions somewhat better than the simple exponentials used by 
Kreslavsky et al. (2000) or Squyres and Veverka (1981) (see section 1.4); we will refer to it as the 
“Hapke-like” form. Then, for example, the reflectance of a rough, faceted or cratered surface, illu- 
minated at incidence angle cos~'(j1,) from its nominal normal vector, within which the unresolved, 
locally flat surface elements are Lambertian, is 1/F' = f(a)Rwo. 

Plane and spherical albedos for rough particles: For a smooth Minnaert surface (S = 0), 
Apr([o) = Ruk-t = (A,(k + 1)/2)u*-! (Equation 1, Table 1), and for this case Apr(u, = 1) = R 
by definition. When S > 0, the value of Apr(j1.) depends on S as well as R. Taking a shadowed 
Minnaert surface as an example: 


R(1 k) 
App (Ho: R, 8) = Ree [ te) yu dude. (7) 


In determining Apr(to, R,S), it is necessary to calculate a for each value of ju, and the integration 
variables (see figure 1). Then, recalling equation 2, A,(R,S) = 2 i Apr(Mo, R, S){odfo. Thus as 
shown in section 1.5, one cannot properly infer the intrinsic surface composition from A, alone - one 
must understand the macroscopic surface roughness in order to back out the “smooth-surface” value 
of R that regolith radiative transfer models or laboratory observations address, which can be related 
to composition. Below, we show how this can be done. 
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Table 1: Relationships between scattering laws and albedos of different types, for macroscopically 
smooth surfaces (S = 0; Cuzzi 1985). When S > 0, the local surface element reflectivity will be 
the product of //F and a function f(a) manifesting surface shadowing, where a is the phase angle. 
The planar flux albedo Ap, and spherical particle albedo A, shown here, and the geometric albedo 
p as well, then decrease with increasing S (figures 5 - 7). We will refer to the values of A, as shown 
above, valid for S=0, as A,,. For completeness we note that the geometric albedo p = A,P(0)/4 (in 
Hapke 1993, p = A,). “Smooth surface” phase functions P,(a@) shown here are normalized in the 
usual way: [ PdQ = 47. For shadowed surfaces, section 1.4 shows that the full phase function is 
simply given by P(a) = f(a)P,(a) (but must be renormalized). Note that the table caption in Cuzzi 
(1985) mislabels R as, essentially, the //F' at normal incidence and reflection, which is true only for 
a smooth Lambert surface. 


1.4 Phase functions 


Traditional scattering techniques (and unresolved observations) treat particles as points, each having 
a phase function P(a) where a is the phase angle. In this situation the incident flux density 7F 
is transformed into an emergent intensity by J(a) = aF'A,P(a). For an unresolved object, the 
observed flux is simply divided by the (known or assumed) solid angle of the object to obtain this 
average intensity. The phase function is independent of azimuth angle about the incident direction, 
and is normalized over solid angle 2 such that { P(a)dQ = 2m J} P(a)sin(a)da = 4m. For isotropic 
scattering, P(a) = 1. 

As a side note, such a treatment of particles as points is acceptable when particles are widely 
separated, which is not a good assumption for most planetary rings. When particles are closely 
spaced and much larger than the wavelength, it is most appropriate to treat photons as scattering from 
individual particle surface elements, each having its own scattering law. Monte Carlo approaches such 
as Salo and Karjalainen (2003), and ray-tracing approaches in general (Porco et al. 2008, Ciarniello 
2015) can handle limb darkening and shadowing of particles in principle, but sometimes the process 1s 
approximated using the “point” phase functions we discuss here. Monte-Carlo or ray-tracing models 
of actual surface facet scattering are not discussed in this paper, but they can be based directly on 
rough-facet results such as presented here. 

Some smooth, flat surface scattering laws can be directly integrated over the particle surface 
to obtain closed-form solutions for the phase function of the particle as a whole. The simplest is 
the Lambert-surface sphere (Salo and Karjalainen 2003): P,(a@) = (8/37)[sina + (7 — a)cosa}. 
A sphere covered by a Lommel-Seeliger surface also has a closed-form phase function: Prs(a) = 
(3/(4(1 — In2)) [1-sin(a/2)tan(a/2)In cot(a/4)]. The phase function advocated by Hapke (1993) 
(equation 10.46, before accounting for the opposition effect near zero phase angle) is essentially an 
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albedo-weighted combination of Pz; and Pg (see section 1.2 for the role of albedo). 
Dones et al. (1993) proposed a power-law phase function of the form: 


P(a) = C,(m — a)” (8) 


where a is the phase angle, n is a positive constant, generally between 2 and 6 (Dones et al. 1993, 
Doyle et al. 1989), C,, is a normalization constant, given by 


2 
Ch = =x (9) 
Dizo Cri 
where 2i+n+2 1 4 
Cri = = & ) (10) 


(2i+n+2)(2i+ 1)! 
We should note that, when converted to magnitudes, the powerlaw phase function does not have 
the near-linear or slightly concave-upwards behavior seen commonly in observed phase curves of 
atmosphereless bodies, while equations (14) and (15) below do. 

Figure 3 shows various phase functions P(q), calculated by numerical integration over the sur- 
face of a sphere covered with elements that locally scatter as the product of a Minnaert law depending 
only on ({Uo, /4), and a shadowing term (equation 6) depending only on a: 


= fla) ED ky, (11) 

To illustrate how easily the roughness effects separate out, we briefly digress on the integration 
as applied to the shadowed case, which is generally applicable to any stipulated //F' function with the 
general form of equation (11). The integral over observable surface elements dA is most easily done 
in photometric coordinates (Minnaert 1961), where ( is photometric latitude and is photometric 
longitude; then dA = cosGd(@d, (on a unit sphere), and: 


a /2 a /2 
P(a) = Ces / | T(dto, H, cos Bd BAX, (12) 
—1/2 J—(n/2—a) 


where C’ is the normalization constant for the smooth-surfaced object. We call out a shadowing- 
related normalization factor cs separately, as it enters only when S > 0, and only in the phase 
function. As noted by Harris (1961), the simplification is that for each integration in photometric 
coordinates at some a, a is constant across the lit face, allowing us to take f(a) out of the integral. 

From solid geometry*, ju, = cos(i) = cosBcos(a — \), and ys = cos(e) = cosGcos,. For the 
Minnaert case this results in 


a /2 a /2 
P(a)= esf(a)c’ f cost 3a3. | cos” \cos*(a — d)dd, (13) 


—7/2 (7/2—a) 


where the shadowing function f(a) comes out of the integral because a is formally independent of 3 
and X (see figure 1), and all other constants are lumped into C’. The integral over @ can be lumped 
with C’ into a new constant C”, because it is also independent of a. The remaining integral over 
can easily be done numerically, or, for the Lambert case k = 1, analytically, providing a handy check 
on the numerics. For the more complicated H-function-derived [/F’, or any other function in general, 


‘The sign convention for \ in Hapke (1993, p. 111) is opposite to ours. 
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S | cste | Cstx | Cstse | CSLSH 
0.0 | 1.00 | 1.00 | 1.00 1.00 
0.2 | 1.23 | 1.16 1.25 1.17 
0.4 | 1.49 | 1.35 | 1.53 1.37 
0.6 | 1.79 | 1.56 | 1.86 1.60 
0.8 | 2.14 | 1.80 | 2.24 1.86 
1.0 | 2.52 | 2.08 | 2.67 2AS 
1.2°| 2.96 | 2.39 | 3.15 2.48 
1.4 | 3.44 | 2.73 | 3.69 2.86 
1.6 | 3.97 | 3.13 | 4.28 3.28 
1.8 | 4.55 | 3.57 | 4.93 ae a 


Table 2: The normalization constant cs for a range of S. Subscript L refers to an underlying Lambert 
phase function (equation 14), and subscript LS refers to an underlying Lommel-Seeliger phase func- 
tion (equation 15). Subscript E refers to the simple exponential f(a) (first expression of equation 6) 
and subscript H refers to the Hapke-like f(a) (second expression of equation 6). 


a tabular lookup can be used to get values of (/1,, 14) uniquely from each (3, A). The constant C” is 
determined by normalizing the smooth surface (S = 0, cs = 1) phase function: [ P,(a)dQ = 4r. 
The fact that f(a) comes out of the integrals makes it possible to write the combined spherical phase 
function for a shadowed Lambert or LS surface (Hapke (1993, equation 12.59); Squyres and Veverka 
(1981)) directly as: 


Ps_1(a) = csf(a)- . [sina + (7 — a)cosa] = cs f(a)P,_1(a), (14) 
or 5 
Ps_ts(a) = csf(a)- Jam" — sins tan 5 Incot = cs f(a) Po_xs(a), (15) 


where the closed-form expressions for P,_; and P,_1s5 are properly normalized in table 1 when S = 
0; when S > 0 a further normalization for each phase function 27cs | f(a)P,-,(a)sinada = 4r is 
needed (see Table 2). So the smooth or rough particle phase functions are determined by the surface 
element scattering law (//F, first row in Table 1 when S=0) for an individual facet. 

Various typical (normalized) phase functions are shown in figure 3, in two formats to make 
the high- and low-phase behavior equally visible. The clusters of solid lines contain curves for Min- 
naert / = 1.0, 1.05, 1.15 respectively; this variation of k produces negligible difference in the phase 
functions. These different clusters are separated by their shadowing parameter S = 0.0 (bottom set), 
0.6, 1.0, and 1.4. A shadowed Lommel-Seeliger particle phase function is shown in the blue dotted 
lines, for S = 0, 0.8 (left), and 1.4 (right). The colored dashed lines are for a Lambert sphere (red), 
and the best-fit powerlaw phase functions cited by Dones et al. (1993) for Europa (cyan) and Callisto 
(Doyle et al. 1989, Dones et al. 1993 (green)). The shadowed-Minnaert functions provide quite a 
good match to observed powerlaw phase functions over the relevant range; the shadowed Lommel- 
Seeliger provides as good a fit at low phase angles, but a less good fit at high phase angles, and lacks 
physical justification for bright particles. The single-scattering shadowed phase function of Lumme 
and Bowell (1981b, equation 27) is of a very similar form to our “Hapke-like” f(a) alone; however 
instead of multiplying a Lambert phase function it is additive. 

For additional reference we show in figure 4 a grid of shadowed-Lambert particle phase func- 
tions (Minnaert k = 1), using either the simple exponential form (solid curves) or Hapke-like form 
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Figure 3: Typical shapes of phase functions for Europa (cyan dashed) and Callisto (green dashed, 
n=3.3 powerlaw) are fairly well matched by different choices of the shadowing parameter S in equa- 
tion 6; Callisto is apparently photometrically rougher than Europa. The red dashed line is an unshad- 
owed Lambert surface. Solid curves are families of Minnaert functions (of which the Lambert case is 
a subset) with different S, and the blue dotted curves are similarly treated LS functions. Left: simple 


exponential f(a) = exp (—Sa); Right: Hapke-like f(a) = exp (—S/tan(a/2)). 


(dashed curves) of f(a) (equation 6), along with best-fitting powerlaw phase functions for Europa, 
Callisto, and Saturn’s ring particles, and phase functions digitized from Pitman et al. (2010) for Rhea 
in two near-IR bands. In their overall shape, the observations seem to favor the Hapke-like f(a). 

It is intriguing that the product of even our simple exponential shadowing function (equation 6) 
and either the Lommel-Seeliger or Minnaert surface scattering law, results in particle phase functions 
that agree with the powerlaw phase functions P,(a@) = cp(a — a)" advocated on purely empirical 
grounds by Dones et al. (1993) for ring particles and other satellite and lunar phase curves, which 
display a concave-upwards dependence on phase angle. Indeed this agreement strikes us as providing 
some degree of physical justification previously lacking for the arbitrarily chosen powerlaw expres- 
sion. Nevertheless, as noted above, our equations (14) and (15) have a shape more closely resembling 
actual phase functions (shape differences are most apparent when converted to magnitudes). 

For completeness, we note in the context of planetary rings, that “particle phase functions” 
backed out of observations of ring / /F to date have invariably assumed the traditional “Chandrasekhar” 
formalism for the scattering behavior of the layer as a whole, or at least a formalism in which P(a) 
factors out (Doyle et al. 1989, Dones et al. 1993, Cuzzi and Estrada 1998, Cuzzi et al. 2002, Bradley 
et al. 2013, Déau 2015, Ciarniello 2015). However, this approach forces P(a) to absorb whatever 
nonclassical behavior is not properly handled by the “Chandrasekhar” formalism. Yn the case of 
Doyle et al. (1989) and Dones et al. (1993), powerlaw phase functions were adopted that automati- 
cally go to zero at high phase angles. However, Déau (2015) obtain P(a) empirically for the B ring 
and C Ring, by backing out the Chandrasekhar layer function. After normalization, these phase func- 
tions are in good agreement with each other; however, both have significant “isotropic” contributions 
and do not go to zero at high phase angles, as do phase functions observed previously for the C Ring 
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black solid: exp shad Lambert, S=0, 0.2, 0.4,... 1.6 

black dashed: Hapke shad Lambert, S=0, 0.2, 0.4,... 1.6 
red and green solid: Europa and Callisto best powerlaws 
blue and magenta solid: Dones et al 1993 A Ring range powerlaws 
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Figure 4: A low-phase-angle closeup of a grid of shadowed-Lambert phase functions using the sim- 
ple exponential form (solid curves) or Hapke-like form (dashed curves) of f(a), for S = 0.0 to 1.4 
in increments of 0.2. We also show best-fitting powerlaw phase functions for Europa (red), Callisto 
(green), and ring particles (magenta and blue) from Dones et al. (1993) and our digitization of ob- 
served phase functions from Pitman et al. (2010) for Rhea (cyan dashed curves) in two near-IR bands: 
1.8m and 2.2,:m (lower), where Rhea is bright in the continuum away from any water ice bands, and 
at 3.6m (upper), deep in a strong ice band where Rhea is very dark. The observations seem to favor 
the Hapke-like f(a). 
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Figure 5: The planar flux albedo App for the Hapke-like f(a) = exp(—S/tan(a/2)) (equation 6). 
The curves show Appr(ji.) determined numerically for shadowed Minnaert surfaces with R = 1 and 
k=1.0, 1.05, and 1.15 (black, red, and green) given four values of S= 0, 0.4, 1.0, and 1.8 (groups from 
top to bottom). The analytical results for S = 0 (see Table 1) are shown by symbols. 


(Cooke et al. 1991, Cooke 1991) or satellite phase functions (Pitman et al. 2010). This might be be- 
cause the “Chandrasekhar” layer scattering formalism used in this analysis had no way to account for 
multiple scattering or high optical depth nonidealities which might be important at moderate to high 
phase angles. Some first order corrections for this nonideality have been presented by Salo and French 
(2010). In the case of planetary rings, future Monte Carlo models which can handle surface element 
scattering by realistic rough-surfaced particles, in realistic closely-packed, “nonclassical” layers, will 
be needed to properly address this distinction. 


1.5 Albedos are a function of roughness 


Planar Flux albedo: In figure 5 we show Ap;-(i,) based on various I/F’ models for Minnaert func- 
tions with R = 1 and k=1.0, 1.05, 1.15 (black, red, green) as shadowed by four values of S = 0.0, 0.4, 
1.0, and 1.8. Recall k = 1 is the Lambert case where, for S = 0, Apr(pi.) = FR. All the curves are 
numerically derived here, but for the top triplet, an unshadowed case or S=0, we also plot the analyt- 
ical solutions (App(u,) = Ru*'; table 1) as symbols. The agreement is quite good. App(ji, = 1) 
decreases with increasing S, even while the surface material is unchanged, because Apr(fi, = 1) 
includes intensity contributions from emergent directions at phase angles away from zero, where 
shadowing increases as S increases. Apr({i.) decreases further as ju, decreases because shadows 
become even more important at smaller .,. While App depends on S, we reserve the parameter 7 to 
refer only to the case S = 0 - a granular but macroscopically smooth surface. 

Geometric albedo: The geometric albedo p = A,P(0)/4 of an entire macroscopic particle 
(see Table 1 and Cuzzi 1985) is insensitive to roughness and shadowing because shadowing is not 
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Figure 6: Correction to the geometric albedo A, for the effect of roughness, as a function of observed 
A, for several Hapke @ values. In determining A,, it is assumed that the opposition surge is due to 
shadow hiding only. Left: opposition surge amplitude 5,=0.5; Right: opposition surge amplitude 
B,=1.0. These models assume a regolith grain single-lobed Henyey-Greenstein phase function with 
asymmetry parameter =-0.3; higher A,, values can be obtained for more backscattering grains (-0.5). 


significant at phase angles less than about 10-20 degrees (Helfenstein and Shepard 2011, Shepard and 
Helfenstein 2011). However, the correction is not entirely negligible, and can be important for very 
bright surfaces, so needs to be addressed if quantitative compositional inferences are to be drawn. 
We show the effect by calculating Hapke’s (1993) smooth surface geometric albedo (Hapke calls 
geometric albedo A,) over a range of regolith grain albedo wy, using his equations 8.22b, 8.25, and 
10.38, and then converting this value of A, to that of a rough, shadowed surface using Hapke’s 
(1993) “U function” (his equations 12.58 and 12.60) for different values of his slope parameter 0. 
We plot the ratio of the rough surface A, to the smooth surface A, in figure 6, as a function of the 
rough surface (observed) A,. Each color represents a different @, and the symbols along each curve 
represent different w,. The corrections are generally small (<20%) except for bright objects with 
very macroscopically rough surfaces. However, Hapke (1993) admits that the theory is questionable 
for very bright surfaces, because illumination of facets by other facets is neglected, and also for very 
large 6 because high order terms have been neglected. So overall, not surprisingly, geometric albedos 
are not sensitive to roughness and shadowing. 

While working with A, (or p) thus has some advantages, the strong opposition effect (espe- 
cially for bright objects such as Saturn’s rings) complicates matters (French et al. 2007, Salo and 
French 2010, Verbiscer et al. 2013; compare figures 6 left and right). There is considerable ongoing 
uncertainty into the nature of the opposition effect(s) and how to properly model it (them); for recent 
discussions see Muinonen et al. (2011) or Déau et al. (2013). Because of this uncertainty and the 
minor influence of the opposition effects on the integrated reflected energy balance (they contribute 
over a very small solid angle), we will emphasize an approach to compositional inference centered 
on the spherical (or single-scattering) albedo A, and avoiding extremely small a. Deemphasizing 
the opposition effect is consistent with our reliance on Hapke (1993) instead of Hapke (2012), which 
differs primarily in adding more complicated treatments of coherent backscattering. 

Spherical albedo: The particle’s spherical albedo A, (equation 2) depends on both its compo- 
sition (as manifested in R, which one can envision measuring or modeling for a grainy but macro- 
scopically smooth surface) and its macroscopic roughness (as manifested in S). This effect can be 
understood by taking a smooth surfaced particle of some composition, and roughening the surface so 
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shadows are seen at all phase angles greater than zero. Except in the opposition surge near a = 0, 
the rough particle appears less bright at all phase angles than a smooth one of the same composition 
because of the shadows, so has a lower value of A, overall. It is not possible to infer compositional in- 
formation from A, alone, without knowing the associated value of S. While Hapke (1993) discusses 
the shadowing effect in chapter 12 of his monograph, he does not provide an explicit correction for 
A, (see section 2 however). 

Figure 7 shows how the particle spherical albedo A, varies with shadowing parameter S for 
both shadowing laws in equation (6). Both shadowed-Minnaert and shadowed-H-function curves 
appear in this diagram. The Minnaert curves all assume R = 1, and the S=0 values would scale 
linearly with 7. In the left panel, we show actual spherical albedos from our model for objects 
with either different Minnaert k, or as computed for different grain albedos w, using the H-function 
approach of Hapke theory (equations 4 and 5). More absorbing material, naturally, leads to lower 
Wg, and thus lower smooth surface (S = 0) spherical albedos; meanwhile, roughness (increasing S) 
decreases these albedos further. 

The curves in the right hand panel show the rough surface albedos normalized by their S = 0 
equivalent value A,,. Surprisingly, they collapse onto a single, essentially albedo-independent family 
for each f(a) which can be approximated by very simple functions €(S) = A;/A; (magenta and 
blue heavy curves in right panel). This behavior makes it very easy to apply the surface roughness 
correction. The ratio £(S) can be thought of as a way of correcting the observed spherical albedo for 
a rough, lumpy particle to the spherical albedo for a particle made of the same granular material if 
its surface were smoothed out. The near-invariance of the €(S) function with grain albedo w, (for a 
given f(qa)) means that the shape of the underlying “smooth-surface” phase function plays a minor 
role, even while it changes from a Lommel-Seeliger law to a Lambert law as w, varies from 0 to 1 
(equation 10.46 of Hapke 1993, and figures 2 and 3). 

All values of observed A, lying on the (closely spaced cluster of) Minnaert curves in figure 7 
(left; dotted lines, for k=1.0, 1.05, 1.15) refer to particles for which a smooth surface would have 
Rk = 1! Yet, depending on roughness, even these highly reflective materials can have much lower 
Apr (figure 5), and A, (figure 7). Uncorrected roughness and shadowing will lead one to overestimate 
the amount of intrinsic absorptive material present. This implication is supported by the H-function 
curves, each for a different grain albedo w, (left panel of figure 7). Lower grain albedos wz, of course, 
produce darker regoliths (lower A,,), and the grain albedo w, is directly linked to particle composi- 
tion. However, A, can also be reduced by increasing S. Note there is no “imposed” normalization 
here of Apr(f1.) or A; associated with S. Only the phase function is renormalized for different values 
of S, as in Hapke (1993). 

The collapse of the curves in figure 7 (right), for all values of intrinsic surface absorptivity 
as controlled by w,, does not imply that any combination of and S is possible. As alluded to 
earlier, if one created a particle from material that truly does not absorb light, both A, and w, must 
be unity (thus also 7), and so it could not have a phase function characterized by strong shadowing 
(a steep powerlaw, or S > 0). If the surface is non-absorbing, photons scattered away from lit surface 
elements would illuminate directly shadowed ones and return to space. As discussed in van de Hulst 
(1957) (section 8.42), the so-called “Schénberg sphere”, a perfectly reflecting, very rough particle, 
has the phase function of a Lambert sphere (that is, k = 1, R = 1, and S =0). 

Some support for this comes from observational, experimental, and theoretical studies. S asso- 
ciated with Europa is lower than with Callisto (figure 3), and the A,-dependent trend of S evident from 
figure 7 and tables 9 and 10 of Pitman et al. (2010) for Rhea (Rhea’s phase curves become steeper 
and more concave upward with decreasing albedo as wavelength is varied). Phase curves for darker 
asteroids are steeper than those for brighter ones (Shevchenko et al. 2016, Penttilaé et al. 2016). Also, 
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Figure 7: Left: Spherical albedo A, as a function of roughness/shadowing parameter S for two differ- 
ent surface shadowing laws f(a): exponential (red; equation 6 left) and Hapke-like (black; equation 
6 right). Dotted curves of each color: Minnaert facet laws using k = 1.0, 1.05, 1.15 and R = 1; solid 
curves of each color: H-function models with w,= 0.3, 0.6, 0.9, 0.99, and 1.0 (wy values increase 
upward). Increased roughness and shadowing decrease the spherical albedo significantly. Right: 
Here, all curves shown in the left panel are normalized by their S = 0 value, defining the function 
A,/Ago = €(S). The curves for a Hapke-like f(a) collapse into an R-and-w,-independent cluster 
(black curves) well represented by £(S) = exp (—S/,/2) (heavy magenta curve); the curves for a 
simple exponential f(a) also collapse (red cluster) and are approximated by £(S) = exp (—S/1.1) 
(heavy blue curve). This difference gives an estimate of the uncertainty in the correction from ob- 
served (rough) to smooth particle albedo. The value of S inferred from observations might differ 
depending on the form of f(a). The fact that the H-function albedo corrections collapse so well 
over a wide range of grain spherical (single-scattering) albedo wz testifies to the insensitivity of the 
technique to the assumed smooth surface scattering law (Lambert or L-S). 
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Figure 8: The scattering asymmetry parameter g (equation 16) for shadowed-Lambert phase functions 
as defined in equation 14 and seen in figures 3 and 4, using the two shadowing functions f(a) defined 
in equation 6. 


Pilorget et al. (2016) see just this effect in the lab (their figure 9 most clearly shows the anticorrelation 
between grain albedo and inferred roughness in the sense of mean facet slope). Finally, Muinonen 
et al. (2011) see nearly Lambertian scattering (S=0) for very rough surfaces, if they are also bright. 
Muinonen et al. (2011) and Pilorget et al. (2016) ascribe their results to multiple inter-grain or facet 
scattering. This behavior is built in to the shadowing theories of Lumme and Bowell (1981a,b) and 
Shkuratov et al. (1999, their equation 28/f) but not that of Hapke (1993, chapter 12). 

We note a crude, but systematic, way to connect arbitrary observational phase functions with the 
rough-smooth albedo correction of figure 7 (right), as long as their scattering asymmetry parameter g 
is known. Figure 8 shows the asymmetry parameter 


Jo. P(@)cosOsinOde 


i 6) — 
aa fv P(©)sinede 


(16) 


(where 0 = 7 — a), as a function of S, for shadowed-Lambert phase functions with the two forms of 
f(a) given in equation 6. For any phase function of arbitrary shape with known g, figure 8 provides 
a value of S which can be used in figure 7 to at least estimate the rough-to-smooth-surface correction 
factor. 

For reference in section 2, we also note at this point the traditional relationship A, = pq (Hapke 
1993, eq. 10.54), where the phase integral q is defined as 


np 
q= 2 f oa sinada. (17) 
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2 Comparison with rough-particle theory of Hapke (1993) 


Hapke (1993) has developed an elaborate, first-principles, essentially geometric-shadowing-based for- 
ward model of the effect of surface roughness on surface brightness. The theory, laid out in chapter 12 
of his monograph, will not be reiterated here and we will simply reference relevant equations. Hapke’s 
development results in multiplicative functions that can be applied to convert “smooth-surface” reflec- 
tivities into “rough-surface” reflectivities. The conversion functions depend on phase angle a (called 
g in Hapke 1993) and a slope parameter 0, and include numerical factors chosen to give good agree- 
ment with a variety of planetary surfaces. One such function K(a,@) converts the phase function 
of a nominally spherical, regolith-covered object with a smooth surface - P,(q) in our notation or 
®(c, wg, 0) in Hapke’s notation? - into the phase function P(a, @) (or ®(a, wy, @)) of an object made 
of the same granular material, but with a roughened surface covered by unresolved facets and craters: 


P(a,0) = P,(a)K(a, 9), (18) 


where Hapke (1993) tabulates (a, @), but says that for phase angles less than 60° it can be approxi- 
mated by . - _ bts we 
K(a,0) = exp [-0.320(tandtan(a/2))/? — 0.526tanOtan(a/2)| (19) 


Note that in the formulation of Hapke (1993) the smooth surface particle phase function ®(a, wz, 0) 
itself is a function of regolith grain albedo w, in a forward-modeling sense, whereas in our model there 
is no explicit dependence. Instead, we work backwards from a phase function that is determined by 
observations, assuming the rough surface shadowing correction operates on a smooth-surface phase 
function that is nearly Lambertian (it could be Minnaert or Lommel-Seeliger with little effect). 

Comparing equations (14) and (18) shows that Hapke’s K(a,@) plays the same role as our 
function f(a) (equation 6): both multiply a smooth surface particle phase function. In our model, 
the smooth particle is assumed to have a composition-independent Lambertian phase function char- 
acterized by a surface element planar flux albedo R and a corresponding spherical albedo A,, (Table 
1). One can then model 7 or A,, separately using traditional Hapke, Shkuratov, or van de Hulst type 
granular regolith radiative transfer models while neglecting explicit roughness effects (eg. equation 
10.51b of Hapke 1993). 

A comparison between our (normalized) results (right panel of figure 7) and those obtained 
from the complete Hapke theory is shown in figure 9. The Hapke results are shown in the symbols, 
one set for each grain albedo w,= 0.1 - 0.99, with 0.99 being the lowest at each value of S. Hapke 
(1993) does not give a direct correction of spherical albedo A, for roughness and shadowing. Thus, 
to calculate a rough, shadowed particle spherical albedo as a function of his roughness parameter 0, 


we start with his smooth particle geometric albedo A, and use its roughness correction function U (6) 
to determine the geometric albedo of a rough particle as described in section 1.5. Both A, and U(@) 
are parameterized by w,. We also determine the phase integral g (equation 17 above) from the rough 
surface phase function, which is (for this test) corrected from Hapke’s smooth surface phase function 


®(a, wy, 0) using his K (a, @) function (eqs. 12.59 and 12.61).° K(a,@) does not depend on wy, but 


SHapke’s ®(a, wy, 0) has a different normalization than our P,(a), being normalized to its zero-phase value. 

©To be more specific, we use Hapke’s equations 8.22b, 8.25, and 10.38 (which assumes 8.56) for A, and quantities 
going into it, all of which enter into his equation 10.47 (with 8.90) for ®(a, wg, 0). We assume an opposition surge due 
to shadow-hiding only, with an angular width h of 0.3, but one could ignore the opposition effect as it contains negligible 
scattered energy. In realistic cases, the grain albedo w, can be a volume average (Hapke 1993, eqn 7.45a) where for a 
given grain with intramixed “pollutants”, refractive indices are obtained using effective medium theory (Cuzzi et al. 2014). 
Hypothetical, nonisotropic regolith grain phase functions (Souchon et al. 2011) can be treated using similarity relations 
(equations 10.25a and b of Hapke 1993), as discused in section 3. None of these complications regarding realistic materials 
are used in this paper, where we simply assume a range of wg, but we note the thread for the convenience of the reader. 
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Figure 9: Comparison of the shadowing effects on particle spherical albedo A,, from our simple 
model (curves, taken from figure 7) and from the forward model of Hapke (1993) (symbols, for wy, 
ranging downwards from 0.1 to 0.99 at each value of S). The approach is explained in section 2. 
Overall the agreement is fairly good; our model is much less sensitive to composition. 


(a, wy, 0) does. We then calculate A, = A,q to determine the spherical albedo of a rough particle. 
In order to make the direct comparison of models shown in figure 9, we slightly distorted the second 
term in the exact form of Hapke’s A’ (a, @) so that its functional form is identical to the first term: 


K(a,@) = exp [—0.326(tandtan(a/2))'/ _ 0.520(tandtan(a/2))'/?| (20) 


Comparing this form with our own f(a) we see that S ~ 0.844V tanJ; this allows us to compare our 
predictions for A,/A,;. (parameterized by S) directly with those we infer from Hapke (parameterized 
by 0). The general overlap between the two solutions is encouraging, given the simplicity of our 
approach and the various approximations in the approach of Hapke (1993). Even the functional 
form of equation (19), which we have distorted in equation (20) to have the same a dependence as 
equation (6), is only claimed as a good fit by Hapke (1993) for phase angles less than 60°. The decent 
agreement seen in figure 9 is probably due to the fact that, for objects with phase functions such as 
these, most of the energy is scattered at low and moderate phase angle. 

Before closing, we note that the effect described here - a necessary correction from rough- 
surface to smooth-surface spherical albedo prior to compositional modeling - may have implications 
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for compositional inferences from spectral band shape modeling as well. Consider an absorption 
band in the observed J/F' of an icy surface, and assume that geometrical and layer effects can be 
properly removed in converting J/F' to A,. Figure 7 shows that surfaces of larger S require a larger 
correction to obtain “smooth-surface” A,,. We have argued that rough particles made of intrisically 
dark, absorbing material may have phase functions characterized by systematically larger values of 
S compared to comparably rough but highly reflective surfaces (end of section 1.5). This suggestion 
is amenable to simple experimental testing. Thus it is likely that the appropriate rough-to-smooth 
correction is systematically larger for more absorbing surfaces (thus, larger in the depths of a band 
than in the nearby continuum). This means that absorption features - even in a single composition 
surface such as the icy surfaces of the rings and moons of the outer solar system - may appear artifi- 
cially deep when observed well away from opposition, causing uncorrected “rough-surface” albedos 
obtained from them to exaggerate the inferred grain size, especially in ice-rich surfaces where the 
spectral contrast is largest. For example, it has been reported that the regolith grain size is larger in 
the A and B rings than in the C Ring (see a recent review by Cuzzi et al. 2009), and even perhaps that 
the A Ring regolith grains are slightly larger than those on B Ring particles. Yet, this is just the sort 
of systematic bias that one might expect from uncorrected shadowing if the A Ring particles were, 
for whatever reason, rougher-surfaced or clumpier. More study is needed to see if the regolith grain 
size effect is real or not. 


3 Cookbook approach to applying the technique: 


(a) In practice, one first obtains the observed (rough particle) spherical albedo A, and phase function 
P(a) by measuring the brightness as a function of phase angle a, from observations over a range of 
geometries (figure 3). For an isolated object, the fit for the phase function is straightforward. For 
a ring system of optical depth 7, this step usually involves fitting the so-called ‘Chandrasekhar’ or 
classical, single-scattering expression 


UP Oca a 


to a collection of observations at different a, 444, and jy values. It should be remembered that any 
deviation of the actual layer reflectivity function from equation 21, including a significant multiple 
scattering contribution or closely packed particles, will lead to a distortion of the inferred phase func- 
tion P(a) (for more discussion see Cuzzi et al. 2002, Salo and French 2010, and section 1.4 above). 

(b) Either way, S is then determined by matching the shape of the particle’s observed phase 
function to equation 14 or 15. If unknown, the phase variation can be assumed by analogy with 
similar objects (introducing uncertainty, of course). 

(c) One next corrects A, for shadowing, using one of the collapsed normalized (blue or magenta) 
curves of figure 7 (right) which we call €(S). The smooth particle spherical albedo As, = As(S = 0) 
is given by A,, = A,/&(S). For the Lambert surface facets underlying figure 7, we can identify R = 
Ago (Table 1). This observationally determined A,, can then be compared with various theoretical 
model values for smooth-surfaced particles of different composition, as described in the following 
steps. 

(d) The first step in modeling A,, is to get the grain albedo w,, which depends only on the grain 
radius and refractive indices; see Hapke (1993, chapter 6), Cuzzi and Estrada (1998, equations 2-5), 
Hendrix and Hansen (2008, equations 2-3), or Shkuratov et al. (1999) as discussed by Poulet et al. 
(2002). The refractive indices can be just taken for the pure material, or an “intramixture” of material 
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within a single grain can be approximated by a volume-fractional blend (Cuzzi and Estrada 1998) or 
obtained more accurately using Effective Mixture Theory (Cuzzi et al. 2014). An average value of 
Wy for an “intimate mixture” of grains having different compositions can be calculated as in Roush 
(1994). 

(ec) The following basic parameters are determined from w, using equations 8.22b and 8.25 of 
Hapke (1993), assuming isotropically scattering regolith grains of albedo w,: 


y=,/l—w, and pea (22) 
fy 

(f) The approaches we use below are based to varying degrees on isotropically scattering re- 
golith grains.’ In general the regolith grain phase function P,(«) is approximated by a Henyey- 
Greenstein phase function: P,(a) = (1 — g”)/(1 + 9? — 2gcosa)*/?, where g is the the scattering 
asymmetry parameter (equation 16 above and equation 6.7 of Hapke 1993). This phase function 
is normalized in the usual way such that for isotropic scattering, P,(a) = 1. The Hapke theory 
does not calculate a grain scattering asymmetry factor g, nor does it include diffraction of the scat- 
tered light in w, as does Mie theory. However, even without diffraction, the reflected and transmitted 
light are already significantly anisotropic for transparent, icy grains and even for the dielectric slabs 
on which the Hapke wz are based (Liou and Hansen 1971, Liou et al. 1983). Poulet et al. (2002) show 
how an effective g can be derived from the model of Shkuratov et al. (1999), or a range of g values can 
simply be assumed (Hapke 1993, chapter 6 and figure 6.9). We model anisotropic grain scattering 
phase functions in several ways (see gl-g3 below). The so-called “similarity relations” convert the 
albedo of a nonisotropic scatterer into that of an isotropic scatterer. That is, the directly calculated 
Wg can be combined with any assumed asymmetry parameter g into an “effective” value w,; for an 
isotropically scattering grain using equation (14) of van de Hulst (1974), as adopted by Hapke (1993, 
his equation 10.25a): 


1— 
uy = 0, (4 ) (23) 
g 
Substituting w,; back into equation 22 above, we get 
l-—w V2 
ee ei 


This “effective isotropic” value of y, and if needed, the associated value of r,, can be used in any of 
the techniques below to get the spherical albedo A,, for a smooth particle covered in regolith grains 
with albedo w, and asymmetry parameter g. The basic similarity relations above (van de Hulst 1974) 
were developed for situations where the diffraction component is included in the particle albedo; 
because this is not true for Hapke’s w,, some restraint must be exercised in allowing g to approach 1, 
because such large values are generally only seen for phase functions that include diffraction. 

(gl) One simple prediction for A,, is based on equation 10.23 or 10.51b of Hapke (1993; they 
are the same). These equations for the spherical albedo of a smooth, regolith-covered particle derive 
from the reflectivity of a semi-infinite layer of isotropic scatterers, and can be written in two slightly 


different ways: 
1- 1-— 
Ao = 722 (1- tn en a eee (25) 
ee 3(1+ 7) 6 


7A more recent model (Hapke 2002) has a more complete treatment of grain scattering asymmetry. The reader is 
invited to explore this approach as a possible refinement. 
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(g2) Cuzzi and Estrada (1998) used a slightly different, semi-empirical similarity technique 
from van de Hulst (1974, 1980): 


(L= 7%) (1 = 0.1394; ) 


Age = ’ 


(26) 


using equation (24) for 7;. We have compared the smooth particle spherical albedos from equations 
(25) and (26), and they agree to better than 1%, not surprising from comparison of their respective 
definitions. 

(g3) A third approach combines models for smooth particle geometric albedo p (Hapke 1993, 
equation 10.38) and smooth particle phase function ®(a,w,,0) = ®,(a) (Hapke 1993, equation 
10.47), using the standard relation A,, = pq (Hapke 1993, equation 10.54) with q given by equa- 
tion (17) above and P(a)/P(0) = ®,(a). The smooth particle phase function &,(a@) combines 
albedo-weighted Lambert and Lommel-Seeliger phase functions to represent the underlying multiple 
scattering component, with an explicit single-scattering term which is g-dependent. The opposition 
effect is neglected (B = 0), because it contributes negligibly to q, but the regolith grain phase function 
P,(a) described above in 3f enters explicitly. Then, 


1 1 
D=Po (; | sre) | 2(P,(0) — 1) (27) 
and 
2 2 
®,(a) = 7 jo ele) 4) ar] E — sin5 tan SIncot > 
2 Ta _ 
| Are — (1 od . (28) 
6p T 


Equation 27 involves approximations in the first term, and equation 28 is somewhat ad hoc (but 
plausible; see table 1 and section 1.4). 

The optimum way to adapt this approach to anisotropic regolith grain scattering is unclear, and 
we have explored three approaches. Probably the most widely used approach (g3a) uses the full 
Henyey-Greenstein regolith grain anisotropic phase function P,(«), for arbitrary asymmetry 
factor g, in the first order term (1 + y)?(P,(a) — 1)/4 above, and makes no further similarity 
adjustment (equation 23) to w,. ®,(a) must be renormalized for each choice of asymmetry parameter 
g; see section 2. As it turns out, this approach leads to results similar to (gl) and (g2) above in 
the range of (g, w,) of most interest to icy bodies. Two other arguable approaches are to make the 
similarity correction to wy, and either (g3b) continue to use the anisotropic single order scattering 
term as above, or (g3c) assume isotropic scattering (g = 0) in the first order term as well. Relative to 
the nominal case (g3a), option (g3b) leads to a stronger dependence of A,, on wy, and option (g3c) 
leads to a weaker dependence. 

The values of A,, predicted by (g1), (g2), and (g3a-c) are shown in figure 10 for the full range 
of w, and g. For moderately high w, and moderately positive g, they are not greatly different (and 
as noted above, gl and g2 are nearly identical). The overall similarity of case g3c to cases (g1,g2) 
is misleading because what really matters to modeling observed spectra is the behavior in the appro- 
priate region of (g, w,) space. For the purpose of this paper we merely note that inferences as to the 
underlying composition of the material making up the regolith grains must be regarded as having a 
model-based uncertainty due to differences between these approaches at the tens of percent level in 
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Figure 10: Five different approaches to calculating A,, as a function of regolith grain albedo w, = 
(0,1) and asymmetry parameter g = (—1, 1), as labeled on the lower left hand plot; see section 3 
for definitions of the approaches denoted g1,2,3a,3b, and 3c. The effect of the anisotropic regolith 
grain phase function appearing explicitly in the first order term (g3a, g3b) is evident in the stronger 
dependence of A,, on g for those cases, especially at low wy. 


Aso over the range (g1,2)-(g3a,b,c). Poulet et al. (2003) also showed that different assumptions re- 
garding the fine-scale distribution of compositional mixtures (“intimate” vs. “intramixed’’) can have 
significant quantitative implications about the average w, and thus relative abundances. 

(h) Following this “cookbook”, one can trace the path connecting observed brightness spectra 
to the underlying composition. Inspection of the equations defining w, (step d, above) shows that 
the regolith grain diameter D and the mass fraction f of intrinsic absorbing material enter only as 
their product fD. An application of this approach to determine particle composition from spectra of 
Saturn’s rings can be found in Cuzzi et al. (2016). 


4 Summary 


Corrections to spherical albedo (also known as single-scattering albedo) of an isolated particle due to 
roughness and shadowing effects are not negligible, based on known phase functions. For an object 
such as Callisto (or Saturn’s ring particles which have comparable phase functions), where S ~ 0.8, 
the albedo correction due to rough surface scattering alone is almost 50%, and the albedo for Saturn’s 
ring particles, for example, without roughness correction, has been estimated at around the same value 
(Dones et al. 1993, Doyle et al. 1989)! There are clearly implications for how much “absorber” one 
can add into the mixture after allowing for roughness, which in turn have important implications for 
the composition of surfaces of various objects. We see the most immediate application of these results 
as to objects that are studied photometrically over a range of phase angles, such as ring particles or 
near-Earth objects (NEOs). 

In this paper we present a very simple way of treating roughness and shadowing effects that is 
based only on actual observations of some object’s phase function. Practical use of the approach is 
outlined in section 3. The results are in decent general agreement with the more complicated forward 
model of Hapke (1993), and contain no explicit dependence on the albedo of the underlying surface. 
The exact functional form of the underlying surface facet shadowing function f(a) makes a little 
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difference (see figures 3 and 5) but once it has been determined (from observations of the phase func- 
tion), it can be directly applied as a correction function for spherical albedo A,; uncertainties in the 
shadowing function are probably smaller than errors induced by neglecting to apply any shadowing 
correction. 
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